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Abstract. 

Cosmic ray acceleration by shocks related with Slipping Interaction Regions (SIRs) in the Galactic Wind is 
considered. SIRs are similar to Solar Wind Corotating Interaction Regions. The spiral structure of our Galaxy 
results in a strong nonuniformity of the Galactic Wind flow and in SIR formation at distances of 50 to 100 kpc. 
SIRs are not corotating with the gas and magnetic field because the angular velocity of the spiral pattern differs 
from that of the Galactic rotation. It is shown that the collective reacceleration of the cosmic ray particles with 
charge Ze in the resulting shock ensemble can explain the observable cosmic ray spectrum beyond the "knee" up 
to energies of the order of 10^^ Z eV. For the reaccelerated particles the Galactic Wind termination shock acts as 
a reflecting boundary. 
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1. Introduction 

The hypothesis that the origin of the cosmic rays (CRs) is 
predominantly the result of diffusive shock acceleration at 
the blast waves of individual supernova remnants (SNRs) 
has the corollary (e.g. Lagage & Cesarskv ll983l Volk 1987] 
Axford [l994fl that this mechanism should only work up to 
particle rigidities (energy per charge) that are a decade be- 
low the so-called knee of the spectrum at rigidities around 
several 10^^ Volts, or somewhat more. The only consis- 
tent way to escape this conclusion is to assume that the 
magnetic field strength at these shocks substantially ex- 
ceeds the typical interstellar values of several fj,G. Such 
high field configurations might be unusually strong stel- 
lar fields in the winds of very massive stars (Berezinsky 
& PtuskinHniHI Volk & Biermann HMSI Biermann [TM5|l . 
or the strong Alfvenic wave turbulence excited at strong 
shock waves, as speculated by Volk l)1984|l and recently 
calculated in a simplified nonlinear model by Lucek & Bell 
iproil and Bell & Lucek (2001). However, also this field 
amplification is limited and for a SN explosion into an es- 
sentially uniform external medium (SNe type la and core 
collapse SNe from stars of zero age main sequence masses 
below about 15 Mq which have little main sequence mass 
loss) knee energies constitute a firm upper bound (Volk 
et al. 12002(1 . Our working hypothesis is therefore that the 
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knee corresponds to the cutoff for the CR sources in the 
disk. Regarding core collapse SN explosions into the winds 
of massive progenitor stars on the other hand, Bell & 
Lucek |PIIT|) have even pondered the possibility for an 
extension of the SNR origin of CRs to energies far above 
the knee, up to 10^^ eV for protons and to 10^^ Z eV for 
heavy nuclei. The observed steepening of the all-particle 
energy spectrum - and presumably also of the energy spec- 
tra of individual nuclei (Kampert et al. I2001|l - above the 
knee at 3 x 10^^ Volt they attributed somewhat vaguely to 
source statistics. Whether such a picture of disk sources, 
with single power law spectra reaching beyond the knee 
up to rigidities of about 10"'^'^ Volt, is consistent with the 
rigidity dependence of the escape from the Galaxy, is an 
open question. In any case no one of the very optimistic 
SNR scenarios aims at an explanation of the all-particle 
spectrum beyond 3 x 10^* eV, the so-called ankle, where 
the spectrum appears to harden again (for a recent sum- 
mary, see Sommers I2001|l . And these ultra-high energy 
CRs are not our concern here either. They may be extra- 
galactic or due to a top-down mechanism of the decay of 
hypothetical superheavy particles left over from the early 
Universe (for a review see Bhattacharjee & Sig H200(7|l . 

Assuming from now on that the CR sources in the 
disk cut off at the knee, we propose here a rather differ- 
ent mechanism to accelerate CRs in the Galaxy to ener- 
gies between the knee and the ankle. This mechanism re- 
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mains nevertheless intimately connected with the mechan- 
ical energy input from stars into the interstellar medium 
in the form of winds and SN explosions. For the process 
to operate we start from the assumption that the Galaxy 
exhibits a significant supersonic mass loss. It is mainly 
driven by the CRs and the hot gas generated in the disk, 
and removes about IMq yr~^ from the Galaxy with a 
terminal speed of several 100 km s^^ that is of the order 
of the escape velocity (Ipavich 1^75, Breitschwerdt et al. 
[TMTiniMI Zirakashvih et al. 1996i Ptuskin et al. [TM7|l . 
Such a Galactic Wind must develop because these CRs 
are tightly coupled to the thermal gas and the magnetic 
field of the interstellar medium and therefore tend to es- 
tablish an essentially infinite scale height above the disk 
due to their ultrarelativistic mean energy. The coupling is 
due to scattering on magnetic fluctuations which the es- 
caping particles excite themselves in the low-density halo 
of the Galaxy. Supersonic speeds are eventually reached 
at distances of about 20 kpc corresponding to the size of 
the disk. 

The new aspect which we emphasize here is based on 
the consideration that even on a large scale this outflow 
must be far from regular in space and time, and that its 
variations should be of the same order as the mean. We 
shall argue that due to the rotation of the Galaxy the 
differences in flow speed will lead to strong internal wind 
compressions, bounded by smooth CR shocks in an ex- 
panding halo gas whose internal energy density is domi- 
nated by that of the CRs from the disk. These CR shocks 
reaccelerate the most energetic particles produced in the 
disk by about 2 orders of magnitude in rigidity. The re- 
acceleration essentially ensures the continuity of the en- 
ergy spectrum at the knee. A fraction of the re-accelerated 
particles will return to the disk, filling a very thick (sev- 
eral tens of kpc) region including the Galactic mid-plane 
rather uniformly and isotropically. Since at these energies 
the propagation in the average wind environment is diffu- 
sive (Ptuskin et al. ll997|l . also their spatial density in the 
disk is about equal to that of the knee particles. We shall 
call the reaccelerated particles here "Wind-CRs", and the 
primary CRs from compact regions in the disk we shall 
call "Disk-CRs", for obvious reasons. These names indi- 
cate the origin of the particles which is not only different 
in time but also in space. 

Galactic reacceleration mechanisms have been fre- 
quently considered in the past. This concerned reaccelera- 
tion of CRs by shocks in the Galactic disk (e.g. Blandford 
& Ostriker TW, Bcrczhko & Krymsky Hl??^ Bykov & 
ToptyginJJ88, Ip & Axford lTMTl Klepach et al. 120001), or 
acceleration of particles at the termination shock of the 
Galactic Wind ( Jokipii & MorfiU T^SBinHZI ) ■ 

Reacceleration in the Galactic disk can not give large 
enough energies of accelerated particles because the size 
and life time of the shocks are quite limited. A different 
case is the acceleration at the Galactic Wind termina- 
tion shock which can be steady state and is very extended 
(about several hundred kpc). The investigation of acceler- 
ation on such a shock results in rather large maximum en- 



forward shock 




Fig. 1. Schematic picture of a Corotation Interaction 
Region in the Solar Wind. The fast stream flows into the 
interaction region through the backward shock, whereas 
the slow stream enters the interaction region through the 
forward shock. Magnetic field and plasma velocity are par- 
allel in the reference frame rotating with the Sun. 

ergies of the accelerated particles of about 10^^ to 10^* eV 
for CR protons (Jokipii & Morfill lT!?^ On the other hand, 
this form of reacceleration faces the difficulty of the ob- 
servation of these particles in the disk. The matter is that 
the condition of efficient acceleration at the termination 
shock coincides with the condition for strong modulation 
of particles inside the Galactic Wind flow. For a diffusion 
coefficient which increases with energy only the highest 
energy particles accelerated at the termination shock can 
be observed near Earth. 

The alternative source of shocks in the Galactic Wind 
flow which we have indicated above are interaction regions 
far above the rotating disk. They are somewhat similar to 
the so-called Corotating Interaction Regions in the Solar 
Wind (Fig. 1), although the differences in their character- 
istics are quite important. Galactic interaction regions are 
due to the fact that much of the high-mass star forma- 
tion in the Galaxy and the associated active regions such 
as superbubbles and OB- associations, including most of 
the SNRs, is concentrated in the spiral arms. These ro- 
tate relative to the gas in a quasi-stationary pattern with 
a related mechanical energy input from stellar winds and 
SN explosions that implies an enhanced wind in terms 
of an increased mass velocity (see also Breitschwerdt et 
al. I2()()2h . In the Parker spiral- type magnetic field struc- 
ture (Parker 119581) that the outffow from a magnetized 
Galactic disk generates (Zirakashvili et al. I199()|l . faster 
wind streams from more active regions will begin to over- 
take slower streams after the wind has formed in the upper 
halo. If the relative flow speed is supersonic, the interac- 
tion will lead to a shock pair like in the Solar Wind. As 
long as the flow speed of the wind is still low for a given ro- 
tation rate of the wind source region, quasi-periodic com- 
pression regions will form that lead only to a series of 
forward shocks after a steepening time, interleaved with 
rarefaction waves. This distinct situation applies to the 
Galaxy, as we shall see. Indeed, the flow time to the sonic 



H.J. Volk, and V.N. Zirakashvili: Cosmic ray acceleration by spiral shocks in the galactic wind 



3 



point, from a reference level ~ 1 kpc above the disk mid- 
plane where the smooth wind flow hydrodynamics is al- 
ready a good approximation, is larger than or comparable 
with the Galactic rotation time. In this sense the Galaxy is 
a fast rotator whereas the Sun is a slow rotator: The wind 
flow time from the base of the corona to the sonic point 
is short compared to the solar rotation time. A second 
difference is given by the fact that the magnetic field lines 
are rooted in the gas of the disk, especially in the massive 
molecular clouds, rather than in its spiral pattern, and are 
therefore only temporarily compressed in the spiral arms. 
The angular velocity of the interaction regions equals the 
spiral pattern angular velocity and the angular velocity of 
frozen-in magnetic field lines equals the Galactic angular 
velocity. 

As a result, the interaction regions propagate as a 
wave phenomenon across the magnetic flux tubes. From 
the point of view of the interaction regions the magnetic 
field lines slip through them. That is why we shall rather 
call these interaction regions Slipping Interaction Regions 
(SIRs). This slip also means that particles that were ac- 
celerated at one of the shocks will be able to leave the 
shock region eventually to escape along a field line, rather 
than only across the field like in the Corotating Interaction 
Regions of the Solar Wind. To this extent the acceleration 
properties of the rotating Galactic halo arc quite differ- 
ent from those of the Corotating Interaction Regions of 
the Solar Wind discussed many years ago by Fisk & Lee 

In this paper we shall consider reacceleration of Disk- 
CRs by SIR shocks in the Galactic Wind flow, the large 
scale propagation of the high energy Wind-CRs in the 
wind-disk regions, and the resulting anisotropics. The pa- 
per is organized as follows: Brief descriptions of CR trans- 
port and of the Galactic Wind flow are given in Sections 
2 and 3, respectively. Numerical calculations of SIR shock 
formation are presented in Section 4. CR reacceleration 
by shocks related with the Galactic spiral structure is in- 
troduced in Section 5 and calculated in Section 6. Section 
7 contains the conclusions. 

2. Cosmic ray propagation in the Galactic Wind 
flow 

A schematic picture of the Galactic Wind flow geometry is 
given in Fig. 2. The azimuthally symmetric flow originates 
in the Galactic disk and extends the frozen-in magnetic 
field. At small heights above the disk the gas velocity is 
perpendicular to the galactic disk. The flow is approx- 
imately radial at large distances from the Galaxy. The 
gas is assumed to be fully ionized in the wind, support- 
ing magnetohydrodynamic waves. These may also be res- 
onantly excited by the anisotropic streaming of the CR 
component. Their amplitudes will ultimately be limited 
by damping through nonlinear wave-particle interactions. 

A complete and self-consistent model of CR propa- 
gation in such a flow has been considered by Ptuskin et 
al. 1(1997,1 . The CR diffusion coefficient is determined by 




Fig. 2. Meridional cross-section of the Galactic Wind flow. 
The direction of the gas velocity is shown by the ar- 
rows. The supersonic gas flow is bounded by the termi- 
nation shock. CR transport is mainly diffusive inside the 
diffusion-advection boundary; outside this boundary it is 
determined by convection in the gas flow (in the dashed 
region). The galactic disk is indicated by the black ellipse. 

Alfven waves generated by the CR streaming instabil- 
ity. CRs produced in the disk - the Disk-CRs ~ diffuse 
along magnetic field lines. This diffusion flux generates 
Alfven waves and they propagate in the Galactic Wind 
flow. The generation of Alfven waves is balanced by non- 
linear Landau damping. The resulting diffusion coefficient 
Dji along the magnetic field does not depend on distance 
and is approximately given by 

Dl ~ 10'^(p/mpc)^^-3cm2s"\ (1) 

Here p is momentum of the particle, rrip is the pro- 
ton mass and 7d is the power law index of the momentum 
distribution of the CR sources in the disk; its numerical 
value is close to 4. The CR spectrum near in and near 
to the Galactic disk is formed as the result of the inter- 
play between diffusion and advection. At small distances 
from Galactic midplane diffusion dominates advection. At 
large distances diffusion is weaker because the gradient 
scale increases and because the magnetic field becomes 
rather azimuthal. Diffusion and advection are comparable 
at some distance i?da • This distance increases with parti- 
cle momentum as p(Td-3)/3 p^^. ^^ic diffusion coefficient 
given by Eq. (1) and a wind velocity u = 300 km s^^, 
^da(lTeV) « 15 kpc, approximately equal to the disk 
radius. 

The supersonic Galactic Wind flow is bounded by a 
termination shock at some distance from the Galaxy. 
This distance depends on the intergalactic pressure which 
is poorly known. Estimates give a value of several hundred 
kpc. One has to expect that the termination shock creates 
strong MHD turbulence downstream, towards intergalac- 
tic space. CR diffusion is strongly reduced in this case and 
should be close to the Bohm limit. For those CRs whose 
propagation characteristics in the Galactic Wind (inside 
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the termination shock) are diffusive, the termination shock 
can therefore be considered as a reflecting boundary. This 
holds up to a maximum particle rigidity that is deter- 
mined by the condition Z?b uR^ where = vrg/3 
is the Bohm diffusion coefficient of particles with gyrora- 
dius Tg and velocity v. Its numerical value is determined 
by the local magnetic field strength. At large distances 
in the then spherically symmetric Galactic Wind flow the 
magnetic field strength is given by 



B « 



Rln{9) sine 
' r[i ' 



(2) 



where 6 denotes Galactic colatitude, Bg is the poloidal 
field strength in the disk, i?g is the Galactic radius and il 
is the angular velocity of Galactic rotation. Then, using 
Eq. (2) , the above estimate Db ^ uRs gives the following 
maximum energy 



lO-iSs-i J \2iiG J 



Ra 



15kpc 



1.2 • lO^^eV, 



(3) 



where Z is the nuclear charge number of the particle. We 
should note that this energy is not small for small 6 be- 
cause the angular velocity of the galactic rotation in- 
creases toward the galactic center. Higher energy particles 
cross the termination shock diffusively. 

2.1. Single power law sources up to -Bmax t^^^ disk? 

Let us assume now that the sources in the Galactic disk 
produce CRs with a single power law spectrum up to 
very high rigidities. If the rigidity of the particles is small 
enough so that the distance to the diffusion-advcction 
boundary Rda{p) < Rs then in the Galactic disk the 
spectrum of such Disk-CRs will be modified relative to 
the source spectrum proportional to P~'^'^R^^{p) ■ This 
happens as long as E < , where E^ is given by 
Rda{Et:) = Rs and corresponds to E^ lO^^Z eV for 
Rs = 300 kpc . Particles with larger rigidities freely dif- 
fuse in the space interior to the termination shock and 
can only escape convectively through this shock in an en- 
ergy independent manner because of the small diffusion 
coefficient beyond. Hence the observed spectrum must be 
the same as the source spectrum. As a result we should 
observe a hardening of the spectrum in the energy range 
Eif < E < ii^max (see Appendix A for an approximate 
analytical calculation of this effect). 

However, this spectral hardening is not observed. We 
therefore conclude that Galactic sources must have a cut- 
off energy smaller than i?* which in turn is essentially 
equal to the knee energy E'kncc — 3 • 10^^ eV. As a conse- 
quence the notion of a single power law source spectrum, 
produced in the Galactic disk, must be ruled out. Its ba- 
sis was the somewhat arbitrary assumption that Galactic 
escape might become even more effective beyond the knee 



(e.g. Hillas 11984(1 . In the light of a consistent physical 
model for the overall CR transport in the Galaxy this 
hypothesis is untenable. We rather believe that reacceler- 
ation in the Galactic Wind inside the termination shock 
should produce particles with energies ii^knoc < E < i^max 
with the same source spectrum as the spectrum observed 
in the disk. In a bottom-up scenario of acceleration from 
lower energies, higher energy particles should mainly come 
from beyond the termination shock. 

2.2. Quantitative definition of Disk-CRs vs. Wind-CRs 

We can now make our qualitative definitions of Disk-CRs 
and Wind-CRs from the Introduction more precise: equat- 
ing E knoc ~ -E* the Disk-CRs are those energetic particles 
with E < -Ekncc, and their sources are in the disk; Wind- 
CRs are the particles with E knee < E < -Bmax whose 
source is of a quasi-diffuse nature in the Galactic Wind. 



3. Spiral structure of our Galaxy, wind shocks 

The spiral pattern of late-type galaxies was explained the- 
oretically by Lin and Shu (1964) as a density wave in the 
distribution of stars in the disk of galaxies. It is believed 
now that our Galaxy possesses a two or four spiral arm 
structure, or even a mixed two-four arm spiral pattern 
(Lepine et al. l2()()l|l . This pattern rotates rigidly with an- 
gular velocity fip. Therefore the dependence of the wave 
amplitude on time t and azimuthal angle ip is given by 
exp(im((p — flpt)), where m is the number of spiral arms. 
The estimates for the angular velocity of the spiral pat- 
tern are very uncertain and controversial. The old value 
is Hp = 13.5 km s^^ kpc~^ (Lin et al. I1969|l . Lepine and 
Amaral H1997|l found fip = 20 km s~^ kpc~^, Lepin et al. 
(2001| argue in favor of Hp — 26 km s^^ kpc^^. On the 
other hand Fernandez et al. ( 2001 ) give the value Hp — 30 
km s~^ kpc~^. In any case, the Sun is not located far from 
the corotation radius at which the Galactic rotation ve- 
locity equals the pattern velocity. The rotation velocity of 
the Galaxy at the Sun's position is approximately Q — 26 
km s^^ kpc^^. 

Since potential CR sources in the disk are primarily 
concentrated in the spiral arms, we can assume that the 
CR pressure in and above the spiral arms is larger than be- 
tween the arms. Hence, the CR-driven Galactic Wind flow 
should be modulated by the spiral structure. The situation 
differs from the Solar Wind case because of the relatively 
fast rotation of the Galaxy. Indeed, the Galactic Wind 
flow time to the (magneto)sonic point is about 100 mil- 
lion years and is comparable with the period of Galactic 
rotation. This means that the modulation by the spiral 
pattern should produce magnetosonic waves, propagat- 
ing in the Galactic Wind flow. At large distance from 
the Galaxy these spiral compression waves propagate in 
the radial direction and approximately perpendicular to 
the wind magnetic field which is by then practically az- 
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imuthal. The radial and azimuthal wavenumbers for these 
waves are 



1 i^s — 1 



(4) 



where tis the radial velocity of the wave. This value is 
approximately the sum of the wind velocity u and the 
phase velocity of the fast magnetosonic waves Cf traveling 
in the radial direction. At large distances we simply have 



Cf 



where 14 and Cg are the Alfven and the sound velocity, 
respectively. 

The nonlinear steepening of magnetosonic waves can 
produce a train of forward shocks at large Galactocentric 
distances. The characteristic distance can be found from 
the following estimate. The inverse steepening time is 
kr8u, where 5u is the velocity perturbation in the wave. 
During this time the wave propagates the distance 5r ~ 
{u + Cf)/{krSu). Using the expression for the radial wave 
number we obtain 



Sr 



[u + ctY 



For u + Cf ~400 km s~^ , 5u ~ 50 km s~^ , and m = 2, 
Sr is about 60 kpc. The velocity perturbation Su is half of 
the velocity jump at the shock formed. We therefore con- 
clude that spiral shock formation is possible at distances 
of 60 to 100 kpc. Numerical results (see below) confirm 
this estimate. These spiral shocks should not be confused 
with the spiral density wave in the Galactic disk. The ra- 
dial dependence of the spiral density wave in the Galactic 
disk will be transformed into a latitude dependence of the 
spiral shocks in the Galactic Wind flow (sec Fig. 2). 

We should underline that the shocks in the Galactic 
Wind have two distinct features in comparison with the 
Solar Wind. 

First of all the shock striicture is stationary in the 
frame of reference corotating with the spiral pattern. 
Material Galactic rotation still exists in this frame except 
at special Galactocentric radii. Therefore these shocks are 
not in corotation with the matter in the Galactic disk but 
rather slip through it and this is why we shall call them 
SIR shocks. Such a particular feature can be very impor- 
tant for diffusive shock acceleration of particles. This can 
be appreciated by comparing the Galaxy with the Sun. It 
is simple to see that the particles accelerated in the shocks 
of Corotating Interaction Regions of the Solar Wind are 
drawn into the Corotating Interaction Region between 
the forward and backward shocks (see Fig.l). In addition, 
energetic particles diffuse predominantly along magnetic 
lines which are immobile in the frame corotating with the 
Sun. Hence, particles that are moving backward to the 
Sun should pass through the shock again, which is difficult 
because of the small diffusion coefficient generated near 
the shock and in the downstream region of the gas flow. 
Therefore it is only possible to observe these particles in 



the space between the Sun and the Interaction Region (e.g. 
in the neighborhood of the Earth) if they diffuse across the 
field which is a very slow process. Galactic SIR shocks, on 
the other hand, are not in corotation with the wind mag- 
netic field lines because these field lines are anchored in 
the Interstellar gas of the Galactic disk. Hence the inner 
parts of field lines become "free" towards the disk (after 
the shock has left them behind and the shock-created tur- 
bulence has subsided). This allows SIR-acccleratcd parti- 
cles to reach the disk along magnetic field lines without 
any necessity to invoke perpendicular diffusion. 

The second feature is the absence of backward SIR 
shocks. Outward moving large-amplitude periodic waves 
steepen into a train of shocks (saw-tooth wave), where 
forward shocks are followed by rarefaction waves, like in a 
gas at rest. 

4. Numerical calculations of SIR shocks 

In this Section we present numerical calculations of SIR 

shock formation in the Galactic Wind flow. We shall use a 
simplified spherical geometry for the Galactic Wind flow 
and investigate the formation of SIR shocks propagating in 
the Galactic equatorial plane. For simplicity we describe 
thermal gas and CRs in a two-fluid approximation with 
two-dimensional magnetohydrodynamic ccjuations includ- 
ing CRs (see e.g. Zirakashvili et al. (1996)): 



dp 

dt 

du 
'dt 



1 

j,2 Qj, 



(5) 



du 
dr 



Uip du 



dip 



r 



9$ 
dr 





1 



dB 



Anp \ r d(p 



(6) 



duip _ du^ u^ du^ u^u 
_ _ _u 



dr r dip r 



dB^ dB,„ u,„ dB,„ „ d u,„ B,„ d 

= —u 



dt 



"^-^+5r^^-^£(™) (8) 



dr r dip dr r r dr 



r dip dr 



dP„ dP„ u^ dP„ 

-w = -^—--— 



(9) 



dr r dip 



(10) 



dPc _ _ My + ^ay dP^ 

dt ^ dr r dip 

-7cPc {r\u + v,)) + + (11) 
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Here p is the gas density, u and B denote the radial com- 
ponents of the gas velocity and the magnetic field respec- 
tively, and are the azimuthal components of the 
gas velocity and the magnetic field, $ is the gravitational 
potential, Pg and Pc are the gas and CR pressures, respec- 
tively, 7g and 7c with 7c < 7g are the adiabatic indices of 
the thermal gas and the CR gas, and Va_ = B/^Airp and 
Va.ip — B^pl yf^Kp are the radial and azimuthal components 
of the Alfven velocity, respectively. 

Eq. Q describes mass conservation. The overall mo- 
mentum balance Eqs. 10 and ( O determine the mass 
motion in the radial and azimuthal directions. Eqs. (|SJ) 
and lO describe the evolution of the frozen-in magnetic 
field. The time evolution of the CR pressure is given by 
Eq. Hll|) . neglecting diffusion in the frame of the waves. 
This is an excellent approximation for the description 
of these large-scale structures since the CR pressure is 
mainly determined by low-energy CR particles produced 
in the Galactic Disk. The transport of these particles is 
completely dominated by advection in the SIR waves. 
Following Zirakashvili et al. ( I19H6|1 we assume that CR 
streaming generates Alfven waves propagating along the 
magnetic field out of the Galaxy. These waves are damped 
and heat the gas. This dissipative effect is described by the 
last term in the equation for the gas pressure, Eq. H10|l . 

The gravitational potential $ was chosen in the follow- 
ing simplified, spherically symmetric form: 

^.*„_£iM2£ + £iM!i,„(i + J_) ,12, 

for r < 100 kpc. Here G=6. 668-10"^ cm'^g~^s~^ denotes 
the gravitational constant. The second term on the r.h.s. 
describes the gravitational potential of the Galactic disk 
and bulge, and the third term is the potential of the 
Dark Matter Halo with radius 100 kpc. We used the values 
Mb.d = 9.97- IQi^'Mo, Mr = 1.07- lO^Afo, i?H = 12 kpc 
from Allen & Santillan l|1991|l . The total mass of the Dark 
Matter Halo amounts to 8.0 • lO^M©. 

The system of equations © to ((TT|) was solved numer- 
ically using an explicit finite-difference scheme. Artificial 
viscosity was included to avoid numerical instabilities. At 
the base level at r = 15 kpc we fixed the gas number den- 
sity Uq = 0.001 cm~'^, the radial magnetic field component 
Bo = 1.0 • 10-6 G, the gas temperature T = 7.25 • 10^ K, 
and the azimuthal velocity u^o = 220 km . Adiabatic 
indices 7g = 5/3 and 7c =4/3 were used. The CR energy 
flux Fc = 7c(m + Wa)-Pc/(7c — 1), modulated by the spiral 
structure, was taken in the form 

Fc^ Fco {1 + 0.5 sin{m{ip^Qpt))) (13) 

which assumes that the power of the Disk-CR sources in 
the spiral arms is a factor of three larger than between 
the arms. For the numerical calculations we used to = 2, 
rip = 30 km s"i kpc^i and Fco = 2.0-10-5 erg cm-^ s"^ 
Numerical results are shown in Fig. 3. One can see 
that SIR shocks with a velocity jump of the order 100 
km s-^ are formed at distances exceeding 50 kpc. These 



SIR shocks form a saw-tooth wave velocity profile. The 
compression ratio CTs of these shocks is about 2.0. In the 
sequel we shall call the entire saw-tooth wave the "SIR 
shock system" . 

We emphasize that our numerical results show the 
presence of forward shocks only. Backward shocks were 
also observed at initial times in the computation. 
Nevertheless the backward shocks were blown out at large 
times when all quantities depend on r and ip — flpt only. A 
similar result was obtained by Cranmer & Owocki H199()|l 
for the formation of Corotating Interaction Regions in O- 
star winds. 

5. Reacceleration of particles by the SIR shock 
system: physics conditions 

We shall consider now the possibilities for standard diffu- 
sive shock wave acceleration in the spiral SIR shock sys- 
tem. A similar consideration was made by Spruit l)1988|l 
for spiral accretion shocks. Due to the fact that the wind is 
primarily driven by the CRs which have a softer equation 
of state and dominate the thermal part of the pressure 
(Zirakashvili et al. 11996(1 . Galactic Wind SIR shocks are 
basically smoothed by the pressure of the "low" energy 
Disk-CRs (compare Drury & Volk ll98l]l : in addition, they 
are almost perpendicular shocks. In this case there is no 
injection of particles from the thermal pool and SIR shocks 
will only reaccelerate Disk-CRs with energies E close to 
^-knce to cncrgics beyond the knee. 

It is interesting to consider the implications of this sit- 
uation. First of all, the SIR shock system will not modify 
the spectral form of the Disk-CRs below the knee. Since 
the average adiabatic compression in the saw-tooth wave 
is zero, there will also be no net adiabatic number or en- 
ergy density increase of the Disk-CRs; the spectrum is 
still determined by the interplay of the sources in the disk 
and the propagation in the wind. Secondly, the energy re- 
quirements for the production of the Wind-CRs by the 
SIR shock system are minimal because no other particles 
are energized in the process. The total escape energy loss 
of the Wind-CRs is given by their convective energy flow 
through a spherical surface at the radial position of the 
termination shock (see below). 

The magnetic field in the Galactic Wind is rather 
strong. As a result, the Mach number of the wind flow 
is not very large. We expect that SIR shocks are not very 
strong either and that the spectrum of the particles ac- 
celerated by a single shock of the SIR shock system will 
be fairly steep. On the other hand, particles with energies 
of the Wind-CRs are diffusively locked inside the termi- 
nation shock and can be continuously reaccelerated by 
the system of SIR shocks. It is well known that in such 
a situation the spectrum of accelerated particles can be 
harder than that due to a single shock (e.g. Blandford & 
Ostriker .1980. Spruit 11^551 Achterberg 1990, Mehose k 
PopeUnH Klepach et al. 20001) • The only caveat here is 
that all those many-shock acceleration theories were based 
on the test particle approximation in a system of shocks 
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Fig. 3. Radial dependencies, taken at one azimuth angle. The values of the radial gas velocity u (thick solid line, in 
units of km s"'^ ) are given on the left abscissa. The right abscissa shows the values of the cosmic ray and gas pressures 
Pc (thing solid Hue) and Pg (dotted line), respectively (in units of 10"^'^ erg cm^'^), the gas number density n (dashed 
line, in units of 10~'^cm~'^), and the total magnetic field tension B"^ /An (dash-dotted line, in units of 10"^'^ erg cm~'^). 
Forward SIR shocks form a saw-tooth velocity profile at large distances in the Galactic Wind flow. 



with a strictly discontinuous velocity profile. We shall also 
follow this approach here, even though the finite width of 
the actual SIR shocks should not only produce a power 
law distribution of Wind-CRs but at the same time pre- 
vent any increase of the Disk-CRs. 

The Parker spiral angle a between the magnetic field 
and the radial direction is given by the expression 

cos a = (14) 
V f^^T"^ sin^ 9 + 

We shall assume that the Galactic Wind flow is 
bounded by a strong termination shock at distance Pg, 
where the wind velocity drops by a factor a, the shock 
compression ratio. Beyond the termination shock the gas 
flow is subsonic and is dynamically dominated by the sum 
of gas, magnetic, and CR pressures. One may assume that 
the Galactic Wind flow is highly turbulent in this down- 
stream region. We shall assume that for all particles of 
Galactic origin, Disk-CRs and Wind-CRs, diffusion is neg- 
ligible in comparison with advection beyond the termina- 
tion shock. 

The scattering of the energetic particles is weaker in- 
terior to the termination shock. Nevertheless, as already 



mentioned in the Introduction, particles with E < Pmax 
can be accelerated at the termination shock if the diffusion 
coefficient is small enough there. They can not be observed 
inside the Galaxy for the same reason. Very high energy 
particles, with E > E^i^^ and a correspondingly high dif- 
fusion coefficient, coming from the termination shock or 
beyond, can be observed in the disk, but can not be ac- 
celerated at the termination shock. 

We proceed now to the calculation of the reacceler- 
ation of Disk-CRs by the SIR shock system with spa- 
tial period L = 2Tr /kr and velocity jump Au. The con- 
dition for efficient diffusive shock wave acceleration is 

^ {us — u)L/Q, where D is the diffusion coefficient in 
the direction of the shock normal. It means that the dif- 
fusion time (Dk^)^ ^imist be large in comparison with 
the advection time (fcr(us — u))~^. In the Bohm limit this 
condition gives a maximum energy of accelerated parti- 
cles that is smaller than the maximum energy given by 
equation However, the SIR shocks are almost per- 
pendicular and therefore the diffusion coefficient in the 
direction of the shock normal can be smaller than the 
Bohm value if the shock-associated turbulence is mod- 
erate. This implies a larger maximum energy than that 
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Fig. 4. Spectral energy distributions (in arbitrary units) of the CR protons in the Galaxy (solid curve) and at the 
termination shock (dashed curve) for the case without reacceleration. 



corresponding to the Bohm limit. A necessary condition 
for this situation to apply is D\\/D^ < v/us (Jokipii 
I1987|l . For our parameters this condition is fulfilled within 
a large dynamic range. However expression © still gives 
the maximum energy for collective reacceleration , because 
at larger energies the particles leave across the terminal 
shock rather than remaining confined in the reaccelera- 
tion region. Perpendicular coUisionless shocks with Mach 
numbers smaller than 2 should have a laminar structure 
(cf. Sagdeev [19661 Forslund & Freiberg fl971fl . Therefore 
a moderate level of MHD turbulence downstream of the 
individual SIR shocks and small CR diffusion coefficients 
in the direction of the shock normals are possible. Drift 
motions in the inhomogeneous regular magnetic field in 
the latitudinal direction can also diminish the maximum 
energy of accelerated particles since the SIR shock sys- 
tem depends on latitude. However, we expect the regular 
magnetic field to be small because the regular poloidal 
magnetic field component is weak in our Galaxy (cf. 
Han & Qiao I1994|l . Galactic Wind streams originating 
from different parts of the Galactic disk should drag out 
random magnetic field and magnetic loops with sizes of 
the order of 1 kpc. These magnetic disturbances become 
strongly elongated at large distances from the Galaxy due 



to the acceleration and spherical expansion of the Galactic 
Wind flow (cf. Zirakashvili et al. 1200 1). Hence, the dis- 
tant wind is filled by almost azimuthal, sign-dependent 
magnetic fields. But Parker's formula H14|l and similar ex- 
pressions for the magnetic field strength remain valid for 
this case as well. CR diffusion in such a field should be 
very anisotropic. Random drift motions and wandering of 
magnetic field lines would produce an anomalous trans- 
port across magnetic lines. 

In the following we shall concentrate on the problems 
of continuity of the spectrum of the reaccelerated parti- 
cles and of the possibility to observe them in the Galaxy. 
For the sake of simplicity we shall neglect anomalous 
perpendicular transport and treat CR diffusion as one- 
dimensional similar to the paper of Fisk and Lee (1980). 
Let us calculate the angle between the magnetic field and 
the shock normal ag. Using and ()14|l we obtain for 
r » u/f2 

/ VLUs \ 

cos as — Tz 1 cos a (15) 

One can easily see that this can be zero for a Galactic 
Wind flow line originating at a particular galactocentric 
radius, which means that the corresponding SIR shocks 
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Fig. 5. Spectral energy distributions (in arbitrary units) of the CR protons in the Galaxy (solid curve) and at the 
termination shock (dashed curve) for the case with reacceleration. 



are purely perpendicular shocks. The acceleration by these 
shocks will be efficient if D|| cos^ as ^ (wg — u)L/6. On 
the other hand, particles accelerated near the termination 
shock can be observed in the Galaxy if D^^ cos^ a > uRs- 
These two inequalities can coexist only in the case 



ilpU 



1 



< 



{ug — u)L n{us — u)i 



6uRs 



imuilpRs 



(16) 



This condition determines the part of the Galactic 
Wind flow filled by effectively accelerating SIR shocks 
which produce CRs that are observable in the Galactic 
disk. They are the Wind-CRs. 

6. Acceleration of the Wind-CRs 

For the sake of simplicity we shall consider CR propaga- 
tion in the Galactic Wind equatorial plane. The evolution 
of the isotropic part of the CR momentum distribution 
function N{r,t,p) inside the termination shock r < Rs is 
given by 



ON 
'dt 



1 9 2n 2 dN 
— — r^Z?!! cos^a— 



dN 
dr 




ln( 



dp' 
V 



N{p') - Njp ) 
\n(p'/p) 



(17) 



Here p is the momentum of the particle, is the cosmic 
ray diffusion coefficient along the magnetic field, a is the 
angle between the magnetic field and the radial direction, 
u is the Galactic Wind velocity, and Q{p) describes the 
source of CR particles at the radial distance ro = 15 kpc 
from the galactic centre. The CR distribution function 
N is normalized as n = 47r / p'^dpN, where n is the CR 
number density. The last term on the r.h.s. of Eq. p7|l de- 
scribes the collective acceleration by multiple shocks (see 
Appendix B for details) . It combines the adiabatic energy 
losses of the particles between the SIR shocks of compres- 
sion ratio CTs with multiple reacceleration at these shock 
fronts (the first and second terms in the round brackets, 
respectively). 

The boundary condition at the termination shock is 
given by the continuity of N and of the fiux density: 



n 2 dN 

D» cos a — — 



^ l\pdN 
a J 3 dp 



0. (18) 
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Fig. 6. Spectral energy distributions (in arbitrary units) ol the CR protons in the Galaxy for different cut-off forms: 
exponential (solid curve), oc exp(— p^/pj^^^^) (dash-dotted curve) and the sharp cut-off (dotted curve). 



Note that the termination shock is a reflecting bound- 
ary in so far as the downstream diffusion coefficient is 
taken to be zero in this equation. 

The acceleration efficiency of the shock ensemble does 
not depend on energy. Since the loss of particles into the 
downstream region of the termination shock is energy in- 
dependent also, we expect a power law spectrum of the 
reaccelerated particles, the Wind-CRs. In order to obtain 
the corresponding spectral index 7 we multiply Eq. (|17|l 
by the volume element Airr^dr and integrate over r from 
zero to Rs- In the steady state, using boundary condition 
H18|l and assuming weak modulation inside the termina- 
tion shock as well as a power law momentum spectrum for 
the Wind-CRs, we obtain for the index 7 the equation: 

Rs(tAu\ i?sCrAu In (1 — 7/7S 



7 1 + 



3Lu 



Lu 



In (Is 



0. 



(19) 



The index 7 is always between 3 and the single shock 
spectral index 73 — 3(Ts /(cTg ~ !)• We obtained numeri- 
cal solutions of Eq. (|17|) for the Galactic Wind velocity 
u=300 km s-\ i?s=300 kpc, cr = 3, m = 2, = 450 km 
s"\ Au = 100 km s-\ f^p = 30 km s"i kpc-\ n ^ 20 
km s^^ kpc^^. We also take 7s = 6 which corresponds 
to the SIR shock compression ratio as — 2. These num- 
bers determine the solution 7 = 4.9 of Eq. ifTi^ . The wave 



amplitude approximately corresponds to the numerically 
obtained value from the last section. If we assume that the 
Galactic disk rotates with a velocity of 200 km s~^ , in- 
dependent of Galactocentric distance, then condition H16|l 
will be satisfied for the Galactic Wind flow that originates 
between 9 and 11 kpc. 

The scattering of the Wind-CRs with energies larger 
than lO^^Z eV on self-generated Alfven waves is ineffec- 
tive, since the mean free path corresponding to the diffu- 
sion coefficient is larger than 300 kpc. Higher energy 
particles can be scattered by elongated magnetic field in- 
homogeneities in the Galactic Wind ffow. The mean free 
path of these latter particles can be of the order of length 
scale of these inhomogeneities. To take this effect into ac- 
count we used the diffusion coefficient 



Dn = Df, 



cX 



cX 

y 



(20) 



For numerical calculations we take the energy independent 
mean free path A = 2r. Also, a spectral index of Disk-CR 
sources 7d = 4.0 and a self-consistent cosmic ray diffu- 
sion coefficient Z?j| — 10^^p/(mpc) cm^s~^ independent 
of r were used. These values approximately correspond to 
those obtained in the self-consistent model of CR propaga- 
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Fig. 7. Calculated differential spectral flux I{E) (in units of ni^^sr ~^s^^GeV^^) of the CR protons (dashed curve), 
helium nuclei (dash-dotted curve), carbon (dotted curve), iron (dash-dot-dotted curve), all-particle (sohd curve) in the 
Galaxy for the exponential cut-off, and the all-particle spectral flux observed by the KASCADE collaboration (empty 
circles). The chemical composition has been fixed at E = 9 ■ 10^'' eV from Fig. 5 of Kampert et al. I|2()()l|l . The data 
for the all-particle spectrum are also taken from Kampert et al. ( I2U(J1|I . 



tion in the Galaxy (Ptuskin et al. 11997(1 . The high energy 
cut-off of Disk-CRs was taken as Pmax = 3 • lO^mpC. 

We used the spectrum of Disk-CRs sources Q{p) ~ 
p"'''"' exp(— p/pniax)- The value of the maximum momen- 
tum of Disk-CRs sources Pmax was taken as Pmax — 
3 • lO^mpC. 

Numerical results are shown in Figs. 4 and 5. Fig. 4 
shows the spectral energy distribution (SED) of CR pro- 
tons at the termination shock (dashed curve) and inside 
the Galaxy (solid curve), respectively, for the case without 
reacceleration. Note that the resulting acceleration at the 
termination shock is weak because of the assumed large 
diffusion coefficient inside the termination shock. Results 
taking reacceleration into account are shown in Fig. 5. The 
high energy power-law tail appears beyond the Disk-CR 
cut-off. The mismatch of the SED by a factor 3 inside the 
Galaxy is explained by the modulation of the accelerated 
particles in the Galactic Wind: the Wind-CRs must dif- 
fusively penetrate into the inner wind region against the 
expanding flow. 



We found that the strength of this mismatch depends 
on the unknown form of the cut-off. This feature is demon- 
strated in Fig. 6 that compares SED of protons in the 
Galaxy for exponential cut-off (solid line), sharp cut-off 
(dotted hue) and for the cut-off oc exp(— p^/p^^^^) (dash- 
dotted line). 

The differential spectral flux calculated for the dif- 
ferent cosmic ray nuclei and the corresponding all- 
particle spectral flux, as well as the experimental all- 
particle spectrum measured by the KASCADE collabo- 
ration (Kampert et al. (|2Q01|), are shown in Fig. 7. The 
KASCADE collaboration has since then presented two ad- 
ditional analyses (Antoni et al. 120021 Roth et al. I200:i|l . 
For the all-particle spectrum they agree within less than 
20 percent with the data shown here. 

In principle we should use the source power of the 
Disk CRs as an observational boundary condition in the 
Galaxy, and then calculate the particle spectrum below 
the knee. Instead, we have fitted the amplitude of the pro- 
ton spectrum in Eq. 1(17(1 to the observations just below 
the knee. This fitted proton source power of 7 • 10^° erg 
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is within 20 % the same as that deduced by Ptuskin 
et aL 0^997). To put this number into perspective, we 
assume an average mechanical energy Egj^ = 10^^ erg 
to be released per SN explosion, and a 10 percent aver- 
age efficiency 8 = 10^^ of conversion of this energy into 
CR energy. The resulting Galactic Supernova rate is then 
vsN « 1/(45 yr). 



7. Conclusion 

This paper can be considered as an extension of our pre- 
vious models for the dynamics of the Galactic Wind and 
for CR propagation in the Galaxy (Breitschwerdt et al. 
[TMnn^SI Zirakashvih et al. 1996, Ptuskin et al.nW7|. 
We believe that now a convincing physical picture for the 
Galactic halo, including gas, magnetic field and CRs ex- 
ists. CR sources in the Galactic disk produce energetic 
particles - the Disk-CRs - which mainly drive the Galactic 
Wind flow with its frozen-in magnetic field. CR streaming 
along magnetic field lines excites Alfven waves. Nonlinear 
damping of these waves results in gas heating up to tem- 
peratures of the order of a million degrees in the galac- 

The balance of 



tic Halo (cf. Zirakashvili et al. Ili 
wave damping and wave excitation determines the level 
of Alfvenic turbulence which in turn determines the scat- 
tering efficiency and hence the diffusion of the CR par- 
ticles. The CR spectrum observed up to the maximum 
particle energies resulting from the source spectra pro- 
duced by acceleration in Supernova Remnants (presum- 
ably E < i?knee ~ 3 • 10^^ Z eV) Can be fully explained in 
this self-consistent nonlinear model. For the explanation of 
the higher energy part of the spectrum we propose a reac- 
celeration mechanism in the Galactic Wind flow for the 
most energetic particles from the disk. The lower energy 
CRs are produced more effectively in the Galactic spi- 
ral arms and therefore should drive a more powerful wind 
above the arms. The interaction of these fast streams with 
the slow streams from the interarm regions can result in 
shock formation at large distances from the Galaxy. These 
shocks are similar to the Corotating Interaction Region 
shocks of the outer Solar Wind. CR particles from the 
disk can be reaccelerated on these shocks by a diffusive 
shock acceleration mechanism to generate the Wind-CRs. 
The discussion of the very simple results shown in Fig. 5 
and Fig. 7 allows us to conclude that the observations do 
not contradict the idea. 

We should underline the very important role of the 
termination shock. The expected strong turbulence asso- 
ciated with it results in an energy independent loss of par- 
ticles with energies smaller than about lO^^Z eV. In prac- 
tice this is the maximum energy for the model considered. 
Higher energy particles require accelerating sources out- 
side the termination shock. Alternatively they may be the 
ubiquitous result of the decay of superheavy particles left 
over from the early Universe. 



Appendix A: Influence of the termination shock 
on high energy particles. 

We consider a spherically symmetric wind flow with ve- 
locity u. Beyond the termination shock at radius Rs the 
diffusion coefficient fulfills the inequality D « uRs if the 
turbulence level is sufficiently high. Interior to this dissi- 
pation layer we have an interplay between diffusion and 
advection. The radius of the diffusion-advection boundary 
i?da(p) is determined by the condition D(j), i?da) W-Rda- 
Inside this zone advection is small and will be neglected. 
The solution of the pure diffusion equation between tq and 
Rda. is 



N{p, r) ^ No+ 

+ (iVda - No) 



dr' 



dr' 



r'^D{p, r') 



(A.l) 



Here TV^ja and Nq are the CR distribution functions at 
the distances i?da and tq (Galactic disk radius), respec- 
tively. The boundary condition at r = ro and the continu- 
ity condition at r = i?da are 



D 



D 



r=ro 



dN 
dr 

and 
dN 

r=R,. ~ 
As a result 



-Q(p) 



p dN 
%~dp' 



Na 



„2 

3^ 



dp' Q{p') 

u J P' RUp') 
p 



Nn=3^ 



r'o fdp' Q{p') 



flda(p) 



u J p' RLip') 

P 



dr' 



r'^D{p, r') 



(A.2) 



(A.3) 



(A.4) 



,(A.5) 



If the diffusion coefficient is inversely proportional to 
(this is the case in the self-consistent model of Ptuskin 
et al. ( 1997|) because of the factor cos^ a), the integral over 
r' in the last expression is determined by the upper limit 
and is comparable with the ffist term. A source momen- 
tum spectrum with index 4.0 and a diffusion coefficient 
proportional to p give the value 14/3 ~ 4.67 for the index 
of the observable spectrum . 

For larger energies the radius Rda{p) is comparable 
with Rs and limited by it. Therefore, for these energies 



OO 



rlQip) 



dr' 



r'^D{p,r') 



(A.6) 



This means that the momentum spectrum becomes 
hard like the source spectrum. 
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If, for example, the diffusion coefficient was indepen- 
dent of r, then the second term can dominate the first term 
for uRs « D « uB^/tq in the last equation. But in this 
case the spectrum inside the Galaxy would be determined 
by diffusion only and one should take for example a value 
of 4.0 for the source index and the diffusion coefficient 
proportional to p^^^'^ in order to reproduce the observable 
spectrum. It becomes again hard for high enough energies 
D » uRl/rf). 

Appendix B: CR acceleration by the SIR shock 
system 

Let us investigate the diffusive acceleration by a periodic 
shock wave. Let the velocity profile be a continuous func- 
tion u{x) for Q < X < L where L is period of the wave. 
Shocks are located at x = 0, a; = L, a; = 2L, etc. At each 
shock the fiow velocity drops from the value ui = u{L — 0) 
to the value U2 = u{0 + 0). 

We shall consider the case of a small diffusion coeffi- 
cient, characterized by the inequality D << uL. In this 
case the particle is convected by the fiow from one shock 
to the neighboring one, where it is diffusively accelerated. 
Subsequently it is convected to the next shock where it is 
accelerated again, etc. 

Let N{p) be the particle distribution function in mo- 
mentum p. It is normalized as / Airp^ dpN (p) = ri, where 
n is the particle number density. Near the shock particles 
are accelerated and N{p) changes correspondingly. The 
relation between upstream and downstream distributions 
is given by diffusive acceleration theory (e.g. Blandford & 
Ostriker lTMiji 







(B.l) 



where 7s = 3CTs/(crs — 1) is the power law index of par- 
ticles accelerated by a single shock with compression ra- 
tio tTs = uxjui. For our purposes it is more convenient 

oc 

to use the Mellin transform F(s) — J dpp^^^N{p). Then 



-Fdown(s) — ^up(s)/(l — s/7s)- Between the shocks a par- 
ticle loses energy adiabatically. The distribution function 

1 /3 

N(p) then transforms into N(pas )■ For the Mellin trans- 
form such a transformation reduces to a multiplication 
with (Ts • Hence, after n acceleration cycles the Mellin 
transform of the distribution function is given by 

F„(s)=Fo(s)ar"/V(l-s/7sr = 

= Fo(s) exp (-n (| In + In (1 - 5/7^)) ) (B.2) 

One can check that this does not depend on the ini- 
tial position of the particle. It can be near the shock or 
somewhere between the shocks. Let us now introduce the 

L 

cycle duration T = J dx/u{x). Then, at the time t = nT, 



the function Fn{s) can be found as the solution of the 
equation 



dFjs) 
dt 



T 



(|lnas+ln(l-s/7s)) 



(B.3) 



with the initial condition F{s)\^^f^ = Fo(s). This equation 
for the Mellin transform F(s) corresponds to the following 
equation for the distribution function N{p): 



dN 
'dt 



1 

T 



Ina f 
3 dp J P' \P 





N{p') ~ N{p) 

Mp'/p) 



(B.4) 



It is simple to verify that the Mellin transform of this 
equation gives Eq. ( IB.3|) . In this procedure the value of 
integral 



dy 

y 



-ay 



-by\ 



ln{b/a) 



a > 0,6 > 



should be used. 

It is important that the adiabatic energy losses de- 
scribed by the first term in square brackets of Eq. ljB.4|l 
enter this equation additively. This derivation is also valid 
in the presence of energy independent additional energy 
losses or escape of particles. The corresponding terms 
should be added to the right-hand side of Eq. (|B.4(I . 

For weak shocks Cs — 1 << 1 the operator on the r.h.s. 
of Eq. (|B.4(I is reduced to diffusion in momentum space. 
The corresponding diffusion coefficient is given by 



_ p^ (Au)^ 
" 18 VsL ' 



(B.5) 



where Vs is the shock velocity and Au — ui ~ U2 is the 
velocity jump at the shock. 

In the case of a saw-tooth wave velocity profile 



U{x) = -U2 + ("1 - "2)-^ 



(B.6) 



T = LlncTg/Au, and the r.h.s. of Eq. ljB.4|l reduces to the 
last term in Eq. (|17|l of the main text. 
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